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A  Statistical  Channel  Model 
for  Adaptive  HF  Communications 
via  a  Severely  Disturbed  Ionosphere 


1.  INTRODUCTION 

Although  HF  skywave  channels  have  been  used  routinely  Tor  several  decades 
and  the  shortcomings  of  this  conventional  communications  mode  have  been  well 
documented,  there  is  now  a  growing  interest  in  developing  "smart  radios"  capable 
of  overcoming  most  of  the  shortcomings  that  now  exist.  Due  to  the  rapid  develop¬ 
ment  of  low-cost  miniaturized  computers,  the  development  of  radios  (with  embed- 
ed  computers),  which  drastically  reduce  channel  outage  rates,  signal  distortion, 
and  interference,  is  continually  becoming  more  feasible  and  attractive.  There 
are  a  wide  variety  of  techniques  that  could  be  employed  by  these  "smart  radios" 
and  these  are  collectively  referred  to  as  Adaptive  HF  technology.  This  paper  will 
categorize  these  techniques  and  outline  a  research  program  based  on  experimental 
measurements,  which  will  assess  the  impact  of  the  dynamic  propagation  charac¬ 
teristics  of  the  ionosphere  on  each  class  of  techniques. 

The  technology  to  adaptively  select  operating  frequencies  that  provide  opti¬ 
mum  received  signal  -to-noise  ratios  (SNR)  has  been  demonstrated  and  continues 
to  become  more  feasible  for  general  use.  The  U.  S.  military  is  currently  employ¬ 
ing  a  technique  based  on  FM/CW  chirpsounders  and  scanning  receivers  (interfer¬ 
ence  monitors)  in  the  AN/TRQ-35  frequency  management  system.  In  a  totally 
separate  effort  Collins  Radio  (Rockwell  International)  has  developed  a  system  that 

(Received  for  publication  2  December  1983) 
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probes  allocated  frequencies  using  an  addressable  FSK  code  to  initiate  contact  and 
monitor  network  connectivity.  Other  notable  technologies  that  have  been  demon¬ 
strated  are  adaptive  equalization  modems,  which  remove  signal  distortion  and 
intersymbol  interference,  and  adaptive  antenna  arrays,  which  reduce  interference 
by  optimal  control  of  the  radiation  patterns.  As  a  first  cut  at  categorizing  the 
techniques  involved,  the  following  defines  those  adaptive  functions  that  hold  the 
most  promise  for  significantly  increasing  the  quality  and  reliability  of  HF  com¬ 
munications  systems: 

a.  Adaptive  Frequency  Management  -  Finding  and  using  the  frequency  that 
provides  the  best  SNR  by  making  real-time  probe  measurements  of  propagation 
conditions  and  automatically  assessing  noise  and  interference  levels. 

b.  Digital  Waveform  Processing  -  This  broad  area  includes  adaptive  narrow  - 
band  interference  (NBI)  cancelling,  spread  spectrum  coding  and  detection  (includ¬ 
ing  frequency  hopping),  anti -distortion  processing  using  adaptive  equalizers,  bit 
interleaving,  and  error  detection  and  correction  coding  of  information  bits. 

c.  Networking  -  The  ability  to  automatically  identify  and  gain  access  to,  via 
a  relay  or  relays,  alternate  geographic  signal  routes  to  bypass  ionospherically 
induced  propagation  outage  conditions  on  the  great  circle  skywave  path.  This  is 
also  called  path  diversity. 

d.  Adaptive  Antennas  -  Phased  array  antennas  capable  of  varying  the  phase 
and  amplitude  characteristics  of  their  elements  to  achieve  attenuation  of  unwanted 
signal  sources  and  enhancement  of  the  desired  signal. 

The  feasibility  of  implementing  these  adaptive  techniques  depends  heavily  on 
the  variations  of  spatial,  frequency,  and  time  correlation  of  the  received  signal, 
having  been  distorted  in  these  dimensions  by  the  propagation  medium.  However, 
the  statistical  distributions  of  several  of  these  critical  parameters  are  insuffi¬ 
ciently  described  in  any  literature  from  past  research.  The  current  program  at 
RADC  (entitled  Adaptive  HF  Propagation  Studies)  is  aimed  at  completing  the 
knowledge  of  these  propagation  dynamics  and  providing  statistical  data  (in  hand¬ 
book  form)  that  will  allow  methodical  optimization  of  adaptive  systems. 

Specifically,  the  RADC  Adaptive  HF  experimental  measurement  program  will 
produce  probability  density  functions  (pdf's)  for  six  critical  HF  channel  param¬ 
eters  and  nine  first  derivatives  of  these  parameters  dealing  with  the  time,  fre¬ 
quency,  and  spatial  response  of  the  channel  under  benign  (quiet)  conditions, 
moderately  disturbed  conditions,  and  severely  disturbed  conditions.  These  six 
parameters  are  shown  in  Figure  1  as  the  Gross  and  Fine  Parameter  for  each 
domain  (time,  frequency,  and  spatial).  The  first  derivatives  are  also  described. 
The  communication  links  to  be  measured  include  several  transauroral  and  trans¬ 
polar  paths  to  insure  that  intense  and  dynamic  conditions  will  be  encountered  on  a 
routine  basis  (see  Figure  2), 


easurement  Objectives.  (See  Section  2  for  complete  desi 


CHANNEL  PROBE/SOUNDING  NETWORK 


Figure  2.  Wideband  Adaptive  Arctic  HF  Testbed 


The  channel  model  being  statistically  characterized  by  these  measurements 
is  a  linear  system  or  "black  box"  model  that  produces  an  output  function  (Y)  for 
any  input  function  (X).  For  these  measurements  X  can  always  be  expressed  as 
an  impulse  function  in  the  time,  frequency,  or  position  (spatial)  domains.  Such 
a  set  of  pdf's  should  enable  specific  trade-offs  to  be  made  for  network  architec¬ 
ture,  transmitter  power,  equalizer  specifications,  type  and  complexity  of  signal 
processing,  waveform  design,  specific  implementation  of  adaptive  antennas, 
channel  probe  techniques,  channel  acquisition  protocols,  etc. 

The  six  characteristic  parameters  and  their  first  derivatives  can  be  organ¬ 
ized  in  a  matrix  as  in  Figure  1  and  will  be  described  in  detail  in  the  next  section. 
These  1-D  parameters  were  chosen  to  minimize  the  ambiguity  inherent  in  multi¬ 
dimensional  functions  and  also  to  make  the  data  applicable  to  the  widest  range  of 
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systems  and  applications  by  allowing  modeling  of  different  combinations  of  system 
parameters. 

2.  MEASUREMENT  OBJECTIVES 

2.1  Time  Domain  Channel  Characterization 

(Subparagraph  numbers  correspond  to  parameter  locations  on  Figure  1) 

2.  l.a  Gross  time  dependence  -  Pdf's*  of  signal  levels  in  dBm  will  be  pre¬ 
pared  from  data  at  0.  5,  0.  8,  and  0.  95  times  the  current  value  of  the  MUF  (maxi¬ 
mum  usable  frequency)  on  each  link.  This  data  will  be  of  some  interest  by  itself 
in  selecting  transmitter  power  and  processing  gain  for  reliable  links,  but  it  will 
also  be  the  basis  for  the  spatial  measurement  of  2.  3.  a. 

2.  l.b  Rate  of  channel  deterioration/ recovery  -  For  a  completely  frequency 
adaptive  system  that  can  track  variations  in  the  MUF  (rates  of  variation  are 
addressed  in  2.  2.b)  the  parameter  of  interest  is  the  rate  of  change  of  path  loss 
on  the  frequency  of  optimum  transmission  (FOT,  for  example,  0.  8  or  0.  95  times 
the  MUF).  Therefore,  these  pdf's*  will  be  displayed  as  dBm/min  for  each  link 
measured,  while  being  automatically  adjusted  (floating)  in  frequency  for  variations 
in  MUF. 

2.  l.c  Fine  time  dependence  -  This  will  be  a  compilation  of  parameters 
relating  to  the  linear  system  impulse  response  of  the  HF  channel  being  meas¬ 
ured.  Data  will  be  displayed  in  3-D  composites,  such  as  Figures  3a  and  3b,  that 
are,  respectively,  variation  with  time,  and  variation  with  frequency  (the  latter  is 
also  known  as  an  ionogram).  Also,  pdf's*  will  be  prepared  showing  the  distribu¬ 
tion  of  response  duration  (typically  0  to  5  msec);  that  is,  the  maximum  time -delay 
spread  spanned  by  the  most  coherent  90  percent  of  the  received  signal  power.  It 
should  be  noted  here  that  the  two  probe  techniques  to  be  implemented  during  this 
measurement  program  are  the  mathematical  equivalent  of  transmitting  a  10-psec 
RF  pulse  envelope  and  measuring  the  pulse  spreading  at  the  receiver. 

They  are  both,  therefore,  a  direct  measurement  of  the  impulse  response 
(also  Green's  function  or  spreading  function)  characterizing  the  system  properties 
of  the  onywave  channel.  Since  there  is  so  much  information  derived  from  these 
measurements,  several  other  parameters  (for  example,  coherence  of  individual 
skywave  modes,  amplitude  relationship  of  identifiable  modes  such  as  IF  and  2F, 
fading  rate  of  individual  modes,  etc.)  will  be  examined.  These  impulse  response 

❖ 

Pdf's  will  be  prepared  from  pre-sorted  data  records  such  that  separate  pdf's  will 
be  available  to  characterize  sunrise,  midday,  sunset,  and  midnight  conditions 
for  each  season,  at  0.5,  0.8,  and  0.95  times  the  current  MUF  and  also  exclusive¬ 
ly  for  benign,  moderately  disturbed,  and  severely  disturbed  ionospheric  condi¬ 
tions. 


measurements  will  establish  the  feasibility  of  building  adaptive  data  rate  modems; 
they  will  also  characterize  attainable  spread  spectrum  processing  gains,  and 
required  channel  equalization  filter  length  (maximum  delay). 

2.  l.d  Rate  of  change  of  impulse  response  -  This  parameter  will  be  obtained 
from  the  Doppler  spread  measurements  available  during  Phase  II  of  our  measure¬ 
ment  program  using  a  Wideband  Channel  Probe.  By  scaling  the  observed  Doppler 
shifts  based  on  the  transmitted  center  frequency  the  rate  of  change  of  the  group 
path  delays  can  be  computed.  Pdf's*  will  be  prepared  showing  delay  length  changes 
in  units  of  usec/sec.  This  data  will  comprise  a  benchmark  for  evaluation  of 
adaptive  equalizer  tracking  (convergence)  rates  and  comparison  of  noncoherent 
modulation  systems  (for  example,  FSK,  DPSK)  vs  coherent  modulation  (for  exam¬ 
ple,  coherent  PSK). 

2.  l.e  Time  delay  vs  frequency  -  The  virtual  reflection  height1  of  HF  sky- 
wave  signals  varies  with  frequency.  Therefore,  broadband  signals  experience  a 
continuous  time  dispersion  and  frequency-hopped  signals  a  discontinuous  shift 
possibly  causing  a  loss  of  synchronization.  Pdf's*  will  be  prepared  showing  delay 
spreads  (of  a  single  propagating  mode)  over  100,  200,  500,  and  1000  kHz. 

2.2  Frequency  Domain  Channel  Characterization 

2.  2.  a  Gross  frequency  dependence  -  Pdf's*  will  be  prepared  showing  the 
observed  values  of  MUF  on  each  link.  This  data  will  determine  optimum  distribu¬ 
tions  of  channel  access  (that  is,  network  control)  frequencies  for  various  network 
protocols  and  geographic  locations.  It  will  also  affect  selection  of  the  VHF  exten¬ 
sion  option  (ability  to  operate  above  30  MHz)  for  various  users  by  indicating  how 
often  and  under  what  conditions  these  frequencies  would  be  useful. 

2.  2.b  Rate  of  change  of  MUF  with  time  -  Pdf's*  will  be  prepared  showing  the 
distribution  of  changes  in  the  MUF  (that  is,  MHz/min).  This  will  indicate  re¬ 
quired  propagationupdate  (that  is,  retransmission  of  probing  or  sounding  signals) 
rates  for  real-time  adaptive  frequency  selection  systems. 

2.2.C  Frequency  response  -  Since  the  various  multipath  components  of  an  HF 
channel  arrive  with  random  phase  relationships  (uniform  distribution  of  phase 
angle)  the  impulse  response  is  complex  (that  is,  must  be  represented  as  I  and  Q 
components).  The  channel's  frequency  response  will  be  measured  directly  during 
Phase  II  by  spectrum  analyzing  the  received  spread  spectrum  (20  or  200-kHz 
bandwidth  pseudo  noise  signal).  The  width  and  depth  of  nulls  in  the  frequency  re¬ 
sponse  will  set  performance  limits  on  narrowband  FSK  performance  and  also  help 
select  parameter  ranges  for  adaptive  equalizers. 


1.  Davies,  K.  ( 1965)  Ionospheric  radio  propagation.  National  Bureau  of  Standards 
Monograph  80,  U.  S.  Dept,  of  Commerce. 
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2.2. d  Frequency  response  vs  time  -  3-D  plots  will  display  variations  of 
2.  2.c  over  periods  of  a  few  seconds  with  a  50-msec  update  rate. 

2.  2.  e  Variations  in  frequency  response  -  The  frequency  spread  of  Doppler 
components  measured  from  a  CW  transmission  will  define  the  rate  of  change  of 
the  complex  frequency  response  of  the  channel.  Since  two  multipath  components 
will  change  from  in-phase  summation  to  out-of-phase  cancellation  for  each  n  rad¬ 
ian  change  in  relative  phase,  the  frequency  difference  between  these  components 
specifies  this  fading  rate  directly  (frequency  is  expressed  in  radians/sec).  There¬ 
fore,  the  maximum  fading  rate  on  a  channel  is  the  inverse  of  the  maximum  Doppler 
spread.  To  maintain  practical  significance  this  data  will  represent  only  the  most 
coherent  90  percent  of  signal  power.  Pdf's*  of  this  parameter  will  directly  spec¬ 
ify  required  channel  equalizer  convergence  rates. 

2.3  Spatial  Domain  Channel  Characterization 

2.  3.  a  Azimuthal  correlation  of  path  loss  -  A  critical  parameter  in  evaluating 
the  effectiveness  of  a  relaying  strategy  for  circumventing  channel  outages  is  the 
degree  of  independence  of  the  various  available  signal  routes.  Variations  in  MUF 
and  path  loss  due  to  range  are  quite  predictable  from  geometry  considerations, 
but  have  a  less  predictable  variation  superposed  over  the  deterministic  variation 
due  to  the  physical  displacement  of  the  ray  path.  This  same  displacement  vari¬ 
ability  (that  is,  the  signal  passes  through  a  different  section  of  the  ionosphere) 
can  be  observed  by  approaching  the  receiver  from  a  different  azimuth  at  approxi¬ 
mately  the  same  range.  This  measures  the  decorrelation  of  ionospheric  proper¬ 
ties  as  a  function  of  distance  between  "reflection  points"  without  the  complication 
of  the  range  variation. 

Pdf's*  will  be  prepared  showing  the  cross  correlation  of  signal  levels  for 
each  of  30°,  60°,  and  90°  azimuth  separations.  The  frequencies  at  which  the 
signal  levels  are  measured  will  again  be  0.  5,  0.  8,  and  0.95  times  the  MUF  (the 
same  raw  data  as  parameter  l.a).  This  normalization  relative  to  the  path  MUF 
will  eliminate  the  effect  of  expected  variations  in  MUF  resulting  from  path  mid¬ 
points  being  in  different  sunlight  illumination  regions.  This  data  base  will 
enable  network  designers  to  realistically  predict  overall  connectivity  for  various 
levels  and  complexity  of  network  architectures. 

2.3. b  Rate  of  change  of  path-to-path  correlation  -  Short-term  (minutes)  and 
long-term  (hours)  variations  in  path  correlation  will  be  investigated  to  determine 
an  optimum  "try  again"  interval  for  occasions  when  no  connectivity  to  the  intended 
receiver  can  be  found. 


* 


See  note  on  page  9. 
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2.3. c  Angle-of-arrival  spread  -  Since  the  ionosphere  is  not  a  smooth  spher¬ 
ical  surface,  but  often  has  extreme  irregularities,  a  skywave  signal  will  contain 
scattered  components  that  traveled  via  a  path  different  from  the  great  circle.  Ray 

trace  simulations  performed  on  computer -generated  model  ionospheres  have  shown 

o  2  3 

arrival  angle  dispersion  of  over  10°.  •  Such  variations  make  conventional  null 
steering  ineffective.  Pdf's*  of  angular  spread  of  the  90  percent  power  scatter 
region  will  be  made  as  well  as  actual  intensity  maps  of  cases  of  particular  interest. 
These  data  will  establish  the  feasibility  of  nulling  the  entire  region  of  the  apparent 
interference  source, 

2.3. d  Change  in  apparent  angle-of-arrival  vs  time  -  For  coherent  scattering 
from  a  region  smaller  than  the  resolution  of  the  phased  array  receiving  antenna 
there  is  only  one  apparent  source  that  could  produce  the  phase  front  observed  at 
the  antenna.  However,  this  source  will  move  rapidly  as  the  many  "  wavelets"  from 
the  scattering  region  sum  to  a  different  phase  front.  If  this  phase  center  can  be 
adaptively  tracked  a  deep  null  can  be  held  on  the  scattering  source.  The  feasibil¬ 
ity  of  building  such  a  system  will  depend  on  the  rate  of  movement  of  that  phase 
center.  This  can  be  measured  by  direction  finding  on  the  source  several  times 
per  second,  using  the  complex  amplitudes  at  each  antenna  element  in  a  phased 
array  to  refine  the  apparent  direction  to  greater  accuracy  than  a  resolution  cell. 

2.3. e  Change  in  apparent  angle-of-arrival  vs  frequency  -  Images  of  the 
scattering  region  power  distribution  will  be  computed  simultaneously  on  two 
"zero  bandwidth"  tones  separated  by  a  variable  frequency  difference.  Cross  power 
correlation  of  these  images  vs  frequency  separation  will  be  investigated.  This 
data  will  establish  the  feasibility  of  nulling  broadband  signals  (for  example,  noise 
jammers). 

2.4  Interference  and  Noiae 

All  the  measurement  objectives  described  above  produce  an  absolute  measure¬ 
ment  rather  than  a  SNR  or  bit  error  rate  probability.  It  was  decided  that  measure¬ 
ments  that  ignore  noise  were  more  universally  applicable  since  unwanted  signals 
depend  greatly  on  receiver  location,  locally  generated  noise,  and  on  the  level  of 
enemy  efforts  to  disrupt  communications  links.  On  top  of  these  manmade  condi¬ 
tions,  the  amount  of  noise  that  propagates  to  a  receiver  depends  on  the  path  loss 

>fc 

See  note  on  page  9. 
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and  MUF  in  the  same  manner  as  the  desired  signal  does.  A  statistical  estimate 
for  this  propagated  component  of  noise  and/or  jamming  can  be  derived  from  our 
propagation  data,  while  local  noise  must  be  measured  at  actual  receiver  sites. 

In  conclusion,  noise  levels  are  critically  important,  but  must  be  treated  on  the 
basis  of  the  specific  system  used  and  its  operating  location  and  mission.  The 
propagation  data  compiled  during  this  measurement  program  will  give  a  good 
estimate  of  signal  levels,  which  when  complemented  by  local  ncise  measurements 
and  an  ECM  threat  model,  will  enable  realistic  design  of  link  margins  and  connec¬ 
tivity  rates. 

3.  MEASUREMENT  APPROACH 

Equipment  fielded  to  make  these  measurements  must  be  highly  reliable, 
easily  operated,  and  require  minimal  attendance,  since  project  personnel  will 
not  be  manning  the  sites.  Also,  since  these  sites  will  be  used  for  on-air  testing 
of  Air  Force  prototype  and  production  systems  the  measurement  equipment  must 
serve  the  dual  purposes  of  gathering  the  research  data  outlined  in  Section  2  and 
performing  real  time  test-support  by  characterizing  the  HF  channels  during 
system  testing. 

3.1  Time  Domain  and  Frequency  Domain  Measurements 

Since  a  linear  system  model  of  the  HF  channel  is  being  employed,  the  desired 
output  of  the  channel  probe  system  is  an  impulse  response  and  frequency  response 
at  a  fixed  center  frequency.  This  fixed  center  frequency  must  also  be  automat¬ 
ically  retuneable  (from  2  to  60  MHz)  to  scan  the  entire  potentially  useful  HF 
spectrum.  Therefore,  two  operating  modes  will  be  employed  using  the  same 
system  and  waveform,  a  sounder  mode  and  a  channel  probe  mode.  The  sounder 
mode  however,  is  just  an  automatically  controlled  sequence  of  channel  probe 
measurements  covering  a  predetermined  range  of  center  frequencies. 

Conventional  HF  channel  measurement  systems  are  the  pulse  sounder  and 
FM/CW  chirpsounder.  As  their  names  imply  they  run  in  a  sounder  mode  but  not 
a  channel  probe  mode.  In  the  sounder  mode  they  only  revisit  a  particular  center 
frequency  once  per  several  minutes;  much  too  slow  to  track  many  of  the  propaga¬ 
tion  parameters  of  interest.  Also,  existing  sounder  systems  only  operate  up  to 
30  MHz.  Both  techniques,  however,  can  be  modified  to  operate  as  channel 
probes. 

3.  l.a  Pulsed  channel  probe  -  This  would  require  the  pulse  sounder  to  stay 
on  a  fixed  frequency  and  transmit  10-  to  lOO-ysec  RF  pulses  at  less  than  100-Hz 
pulse  repetition  frequency  (PRF).  The  pulse  length  determines  the  time  resolution 
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of  the  measurement  (100  psec  is  sufficient  to  evaluate  a  conventional  SSB  channel), 
while  the  low  PRF  is  required  to  eliminate  pulse-to-pulse  ambiguity  from  multi- 
path  returns.  A  100-Hz  PRF  provides  10  msec  of  unambiguous  impulse  response. 
A  complex  amplitude  response  can  be  measured  by  processing  an  I  and  Q  channel. 
Due  to  channel  induced  Doppler,  coherent  averaging  in  excess  of  100  msec  is  not 
feasible.  Intra-pulse  phase  coding  allows  more  processing  gain,  but  this  tech¬ 
nique  is  the  basis  for  the  Direct  Sequence  (DS)  and  matched  filter  systems  covered 
in  paragraphs  3.  1.  c  and  3.  1.  d. 

3.  l.b  FM/CW  channel  probe  -  The  "chirp  channel  probe"  typically  sends  a 
train  of  pulses  (up  to  100  percent  duty  cycle),  each  of  which  is  swept  in  frequency 
through  the  fixed  frequency  of  interest.  The  total  chirp  bandwidth  is  20  to 
100  kHz,  while  the  PRF  is  again  less  than  100  Hz.  Since  the  chirp  waveform  is  a 
type  of  intrapulse  coding  the  receiver  provides  processing  gain,  which  reduces 
the  requirement  for  transmitter  power  (for  example,  100  Watts  is  sufficient  for 
most  forward  scatter  applications).  Demodulation  at  the  receiver  is  performed 
by  a  chirping  local  oscillator  (LO),  which  sweeps  in  frequency  at  precisely  the 
same  rate  as  the  transmitter.  Therefore,  all  Intermediate  Frequencies  (IF)  are 
fixed  frequency  signals  with  a  spread  of  spectral  components  offset  in  frequency 
in  proportion  to  the  propagation  delay  experienced  by  various  ray  components  in 
the  received  signal.  The  noncoherent  channel  impulse  response  (that  is,  magni¬ 
tude  of  the  complex  response)  can  be  represented  by  the  power  spectrum  of  this 
IF  signal.  A  conventional  HF  chirpsounder  as  will  be  used  in  Phase  I  of  the  RADC 
program  does  exactly  the  same  thing  except  that  the  bandwidth  covered  is  28  MHz 
rather  than  20  to  100  KHz,  and  the  PRF  is  1/5  min. 

3.  l.c  DS  channel  probe  -  The  next  alternative  to  be  considered  is  a  synchro¬ 
nous  DS  spread  spectrum  probe.  All  DS  systems  require  a  known  spreading  code 
p(t),  preferably  with  very  low  autocorrelation  sidelobes,  that  is,  ideally  this  code 
would  have  the  following  property: 


p(t)  p(t  -  d)T  =  J  pit)  p(t  -  d)  dt  = 


0  for  d  *  0 


1  for  d  =  0 


where  the  integration  period,  T,  indicates  an  ave  age  over  an  interval  T  =  one 
code  length.  For  very  long  codes  small  partial  correlation  responses  (that  is, 
time  sidelobes)  become  insignificant  since  they  will  probably  fall  below  the 
receiver  noise  floor,  making  a  pseudo-random  DS  code  (which  is  not  optimized 
for  time  sidelobes)  acceptable.  The  pseudo-random  codes  are  preferable  for 
secure  communication  systems  since  there  are  nearly  an  infinite  number  of 


practical  code  combinations,  but  for  a  propagation  probe  system  we  can  select  a 
unique  code  with  optimal  correlation  properties.  One  such  code  is  a  periodic 
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maximal  length  sequence  of  2  -  1  bits  per  code  repetition.  This  sequence  is 
trivial  to  generate  using  a  TTL  shift  register  and  one  exclusive  OR  gate  as  shown 
in  Figure  4. 4  5  The  code  rate  and  thus  its  bandwidth  are  variable  simply  by  chang¬ 
ing  the  clock  frequency  fc. 


Figure  4.  Maximal  Length  Code  Generator 

5 

The  length  of  the  code  depends  upon  which  shift  delays  are  chosen  to  feed¬ 
back  to  the  input  giving  a  code  length  N  =  2n  -  1,  where  n  is  the  most  delayed 
shift  involved  in  the  feedback  path.  (Note  that  T  =  N/f,. )  The  auto -correlation  of 
this  sequence  is 

T  -1/N  for  ;j  *  k 

p(j)  p(k)T  =  J  p(j)  p(k)  = 

0  1  for  j  =  k  .  (2) 

For  instance,  if  n  =  10  (N  =  1023)  the  time  sidelobes  of  this  code  are  over  60  dB 
below  the  peak.  This  code  can  be  transmitted  by  bi-phase  modulating  a  carrier, 
(PSK  or  phase  reversal  modulation)  producing  s(t),  at  the  center  frequency  to  be 
probed  resulting  in  the  transmitted  spectrum  of  Figure  5.  The  transmitter  sys¬ 
tem  can  be  quite  simple  as  illustrated  in  Figure  6a  (our  actual  transmitter  system 
is  shown  in  Figure  7,  the  receiver  system  in  Figure  8).  The  received  signal  is  a 
delayed  replica  of  the  transmitted  signal  except,  in  the  case  of  multipath  channels, 
for  the  addition  of  superposed  replicas  having  propagation  delays  (tj)  and  various 

4.  Brookner,  E.  (1977)  Radar  Technology,  Artech  House,  Inc.,  Chapter  8. 

5.  Butler,  F.  (1975)  Pseudo  random  binary  sequence  generators,  Wireless  World, 
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amplitudes  (a^.  Therefore,  using  a  linear  system  channel  model  the  general 
expression  for  the  received  signal  is: 


r(t)  =  /  Channel  (u)  s(u  '  0  du 

=  J  Z  ai  «(«  ■  Ti>  ®(u  ■  ^  du 

i 

=  £  ai  ■  Tj)  sin  [uct  +  $.(t)l  . 

i 


Figure  5.  Spread-Spectrum  Waveform  Spectrum 
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Figure  7  Channel  Probe  Receiver  Site 


Figure  8.  Wideband  Channel  Probe  Transmitter 


This  expression  reflects  the  fact  that  on  a  multipath  channel  the  carrier  is  a 
summation  of  sinusoids  of  random  phase,  the  received  envelope  may  be  a  summa¬ 
tion  of  delayed  code  modulation  envelopes,  and  each  multipath  component  may 
have  an  arbitrary  amplitude.  Since  the  multipath  phases  are  random  a  quadrature 
demodulation  technique  is  required  and,  saying  this  much,  we  will  eliminate  the 
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carrier  term  from  future  expressions.  Due  to  the  correlation  properties  of  the 
code,  the  output  of  the  demodulator  and  low -pass  filter  in  Figure  6a  is  vanishingly 
small  unless  a  controllable  synchronization  delay  (A)  is  set  to  match  one  of  the 
propagation  delays  on  the  channel,  that  is. 


A2  =  ai  p(t  ‘  Ti)  P  (t  ‘ 


a.  for  Tj  =  A 


0  for  r-  *  A 


Therefore,  by  successive  measurements  at  all  delay  settings  of  interest 
^channel^  =  A2^  the  imPu^se  response  of  the  channel  can  be  measured  as  the 
demodulated  output,  hi^)  (a  DC  voltage  level)  for  each  time  delay  Aj.  The  major 
disadvantage  of  this  technique  is  that  there  is  no  assurance  that  the  channel 


response  measured  at  time  A^  will  still  be  valid  when  A^  is  being  measured. 
Since  the  rate  of  change  of  the  impulse  response  parameters  is  one  of  the  meas¬ 


urements  required  of  the  system,  such  an  uncertainty  cannot  be  tolerated.  The 
following  technique  has  been  selected  to  remedy  this  and,  as  shown  later,  makes 
all  delay  measurements  simultaneously,  but  is  the  mathematical  equivalent  to 
synchronous  DS  detection. 


3.2  Matched  Filtered  Detected  Spread  Spectrum  Probe 

The  transmitted  signal  will  be  identical  to  that  of  the  previous  technique 
except  that  the  code  will  be  modulated  onto  the  carrier  differentially  (DPSK) 


rather  than  directly,  allowing  the  received  signal  to  be  nonsynchronously  phase 
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detected.  Nonsynchronous  detection  eliminates  the  effect  of  Doppler  shifting  on 


the  phase  of  the  baseband  output.  Therefore,  the  baseband  can  be  coherently 
averaged  for  on  the  order  of  1  min,  where  synchronous  detection  would  be  limited 
to  about  100-msec  averaging  duration.  The  DPSK  demodulator  will  produce  (for 
a  multipath  channel)  a  baseband  output, 


b(t)  =  £  a.  p(t  -  Tj) 


where  a^  is  the  amplitude  of  the  channel  response  at  time  delay  Tj.  This  output 
will  then  be  sampled  into  a  digital  waveform  analyzer  where  it  will  be  coherently 
averaged  as  many  times  as  required  to  provide  the  desired  SNR.  In  order  to 
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determine  the  channel  impulse  response,  the  waveform  analyzer  can  be  directed 
via  its  interface  bus  to  perform  a  convolution  of  this  sampled  data  with  a  digital 
matched  filter.  Note  that  the  convolution  can  be  performed  by  complex  multipli¬ 
cation  of  the  Fourier  transformed  array  by  the  transform  of  the  code  or  it  can  be 
performed  by  direct  digital  convolution.  The  matched  filter  has  an  impulse 

a 

response  that  is  the  maximal  length  sequence  code  reversed  in  time.  That  is. 

hm(t)  =  •  (5) 

The  waveform  analyzer  will  perform  a  digital  convolution  of  this  matched  filter 
with  the  sampled  receiver  baseband  resulting  in  a  matched  filter  output. 


f 

i 


ym(t)  =  J  b(t)  hm(t  -  u)  du  =  J  Y.  ai  P<u  ■  Tj)  hm(t  -  u)  du  . 


s 


Using  Eqs.  (5)  and  (2)  we  see  that  y(t)  is  itself  a  correct  representation  of  the  non¬ 
coherent  impulse  response  of  the  skywave  channel  (l.c  in  Figure  1), 


y(t)  =  hchannel(t)  =  I  ai  /  P<u  -  Tj)  p(u  -  t)  du 


=  I  a.  6<t  -  Tj)  . 


Therefore,  simply  by  waiting  until  t  =  Tj  the  output  of  the  matched  filter  will  auto¬ 
matically  (that  is,  no  synchronization  or  signal  present  indication  is  required)  be 
aj.  The  moment  that  the  end  of  the  last  code  bit  is  transmitted  is  defined  as 
t  =  0.  Since  a  periodic  code  is  being  used  it  is  important  that  there  is  only  one 
t  =  0  during  the  entire  duration  of  the  channel  response  such  that  all  a.'s  and  tAs 
can  be  referenced  unambiguously  to  the  proper  t  =  o.  Therefore  both  of  our 
operating  modes  have  code  lengths  in  excess  of  10  msec. 
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A  computer  simulation  of  a  coherently  demodulated  PSK  multipath  RF  signal, 
which  is  matched  filter  processed  by  a  127  length  (N  =  2^  -  1)  maximal  length 
code,  is  shown  as  Figure  9.  The  RF  signal  was  time  dispersed  by  discrete 
multipath  delays  and  a  0-dB  SNR  was  created  by  adding  a  random  number  to  each 
sampled  value.  This  simulation  demonstrates  the  tremendous  processing  gain 
offered  by  the  matched  filter,  and  also  demonstrates  the  effect  of  random  phase 
components  on  a  synchronously  detected  signal. 


INPUT  SNR'OdB  l*t  MODE  vt  LO  PHASE* 0° 
CODE  LENGTH -127  2nd  MODE  v*  LO  PHASE*90° 


MAGNITUDE  OF  IMPULSE  RESPONSE  (I2 +  QZ) 1/2 


0  -  -A  ■  «  «  — -  ,  -  .  -  -  I2.7m»«c 

MAGNITUDE  RESPONSE  IN  QUADRATURE  CHANNEL 

Figure  9.  Spread  Spectrum  Processing  Gain 

The  Fourier  transform  of  the  sampled  IF  signal  is  also  of  interest  in  the 
planned  measurement  program  since  it  directly  measures  the  channel  transfer 
function  (2.  c  in  Figure  1).  Sequential  measurements  of  the  IF  spectrum  (for 
example,  each  50  msec)  will  characterize  the  rate  of  change  of  the  channel  trans¬ 
fer  function  (2.  d  in  Figure  1),  since  the  spectrum  of  a  system  output  with  white 
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noise  input  is  the  transfer  function.  The  dynamic  nature  of  these  selective  fades 
in  the  transfer  function  is  shown  in  Figure  10  in  a  4-sec  sequence  of  a  phase- 
modulated  SSB  data  modem  baseband  (same  signal  as  Figure  5), 

The  availability  of  both  time  domain  and  frequency  domain  data  from  the  same 
transmitted  waveform  and  the  simultaneous  processing  of  all  time  delays  in  the 
impulse  response  are  the  major  advantages  of  matched  filter  processing  over  the 
other  techniques.  The  fine  spatial  measurements  required  for  objectives  3.c, 

3.d,  and  3.e  will  be  accomplished  at  a  later  time  under  Phase  3  of  the  RADC 
Adaptive  HF  program.  Therefore,  a  specific  description  of  measurement  tech¬ 
niques  will  be  described  in  later  reports.  The  spatial  parameters  were  included 
in  Figure  1  for  completeness,  but  since  the  hardware  is  still  in  trade-off  analysis 
it  would  be  premature  to  try  to  report  on  findings. 


CHANNEL  TRANSFER  FUNCTION  ASCENSION  IS  -RATRtCK  AFB.  23.5 MHz 
4  CONSECUTIVE  SECOND,  IOIO  EST,  I7-FEB-82 


Figure  10.  Variability  of  Channel  Response 


4.  PROGRAM  PLANS 


The  program  objectives  outlined  in  Figure  1  are  all  to  be  accomplished  during 
three  program  phases.  The  approximate  operational  periods  for  these  phases  are 
Phase  I.  1983-1985;  Phase  II,  1984-1987;  Phase  III,  1984-1987  with  Phases  II  and 
III  possibly  continued  or  reactivated  in  1989-91  to  characterize  a  more  active  sun¬ 
spot  period. 


4.1  Phme  I  —  Azimuth  Diversity  Experiment 

Since  RADC's  propagation  branch  owned  or  had  access  to  three  chirpsounder 
transmitters  and  one  receiver.  Phase  I  was  initiated  during  1982  to  collect  what 
data  could  be  measured  with  these  systems.  The  measurement  objectives  fully  or 
partially  met  by  the  chirpsounder  network  are,  from  Figure  1:  1.  a  (only  below 
30  MHz),  l.b  and  l.c  (below  30  MHz),  1,  d  (no  short-term  tracking,  only  capable 
of  15-min  updates),  l.e,  2.  a,  and  2,  b  (15-min  updates),  3.  a  (no  data  on  true 
one-hop  channels),  and  3.b  (15-min  updates). 

Sounder  transmitters  have  been  installed  at  Ava,  N.  Y. ,  Grand  Forks,  N.  Dak. 
and  Barter  Island,  Ark.  The  transmitters  operate  sequentially  for  5  min  each  in 
a  15-min  cycle.  All  signals  are  received  and  recorded  at  Thule  AB,  Greenland. 
The  sounders  transmit  a  2-  to  30-MHz  signal  swept  in  frequency  at  100  kHz/sec. 
The  demodulated  baseband  is  recorded  on  analog  tape  and  spectrum  analyzed 


during  processing  to  produce  various  amplitude  spectral  lines  representing  the 
propagation  delayed  components  of  the  received  signal  [E  a^  9(t  -  r^)  from 
Section  3 1 . 

4.2  Phase  II  —  Wideband  Channel  Probe  Network 

DPSK  spread  spectrum  signals  will  be  transmitted  at  the  500-Watt  level  at 
the  sites  marked  T  in  Figure  2.  Three  modes  will  be  operated: 

Mode  1  -  10-kbps  length  127  binary  pseudo-noise  code 
Mode  2  -  100-kbps  length  1023  binary  pseudo-noise  code 
Mode  3  -  CW  tone  to  measure  Doppler  spread. 

The  mode,  center  frequency,  dwell  time,  signal  processing  steps,  data  for¬ 
mat,  and  process  sequence  will  be  controlled  by  a  programmable  instrument 
controller.  Therefore,  the  operating  frequencies,  averaging  intervals,  Doppler 
resolution,  sequencing,  etc.,  can  all  be  changed  at  will  by  calling  the  site  and 
uploading  a  new  program.  The  data  collected  at  the  receiver  sites,  marked  R 
in  Figure  2,  will:  (a)  satisfy  all  the  objectives  identified  in  Phase  I  above; 

(b)  eliminate  the  measurement  limitations  noted  in  parentheses  in  Phase  I,  and 

(c)  fulfill  objectives  2.  c,  2.d,  and  2,e. 

4.3  Phase  III  —  Spatially  Adaptive  Antenna  Project  (SAAP) 

The  current  plan  for  this  phase  of  the  Adaptive  HF  effort  is  to  field  a  single 
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large  array  (approx  1  km  )  of  64  elements  capable  of  making  complex  amplitude 
measurements  on  each  element  at  four  center  frequencies.  Signal  sources  will 
be  four  transmitters  all  within  one-hop  skywave  range  of  the  array,  which  will  be 
capable  of  transmitting  two  "zero  bandwidth"  tones  separated  in  frequency  by  a 
variable  amount.  From  the  complex  amplitude  measurements  on  the  array  the 
source  field  (that  is,  secondary  scattering  sources  at  the  ionospheric  "reflection" 
point)  will  be  computed  and  thereby  the  degree  of  spatial  spreading  can  be 
measured.  Since  the  array  is  2-D  there  will  be  resolution  in  elevation  as  well 
as  in  azimuth.  This  system  will  enable  measurements  of  3.c,  3.d,  and  3.e. 

5.  CONCLUSION 

Most  of  the  techniques  included  under  the  general  heading,  Adaptive  HF,  are 
readily  implemented  on  line-of-sight  or  wireline  channels,  however  the  HF  chan¬ 
nel  presents  several  new  technical  challenges  because  of  the  dynamic  properties 
of  the  propagation  medium.  In  order  to  select  feasible  technologies  to  incorporate 
into  modern  military  radios,  the  magnitude  and  rate  of  change  of  these  dynamic 


properties  must  be  thoroughly  characterized.  Since  there  have  been  no  high- 
resolution  short -interval  measurements  of  complex  impulse  response  made  under 
disturbed  ionospheric  conditions  there  is  a  need  for  experimental  data  addressing 
the  time  domain,  frequency  domain,  and  spatial  domain  response  of  the  iono¬ 
spheric  skyway  channel.  The  experimental  program  outlined  in  this  paper  prom¬ 
ises  to  fill  this  gap  in  our  knowledge  in  a  way  that  can  be  readily  applied  to  com¬ 
munication  systems  analysis  and  design. 


8.  Thomas,  John  B,  (1969)  An  Introduction  to  Statistical  Communication  Theory, 
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Acronyms 


decibels  relative  to  a  milliwatt 

direct  sequence  (code  modulation  technique) 

frequency  optimum  transmission 

intermediate  frequency 

local  oscillator 

maximum  useable  frequency 

probability  density  function 

pulse  repetition  frequency 

radio  frequency- 

signal  to  noise  ratio 


MISSION 

of 

Rome  Air  Development  Center 

PAVC  plan&  and  execute*  research,  development,  test  and 
detected  acquisition  programs  in  support  of  Command,  Control 
Communication*  and  Intelligence  (C3 1)  activities .  Technical 
and  engineering  support  within  area*  of  technical  competence 
l*  provided  to  BSD  Program  Offices  ( PO *)  and  other  ESP 
element*.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  control,  sur¬ 
veillance  o f  ground  and  aerospace  objects.  Intelligence  data 
collection  and  handling.  Information  system  technology, 
Ionospheric  propagation,  solid  state  science *,  microwave 
physic*  and  electronic  reliability,  maintainability  and 
compatibility . 
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